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ABSTRACT 
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S-H ' Context. Ammonia is one of the best tracers of cold dense cores. It is also a minor constituent of interstellar ices and, 

as such, one of the important nitrogen reservoirs in the protosolar nebula, together with the gas phase nitrogen, in 
the form of N2 and N. An important diagnostic of the various nitrogen sources and reservoirs of nitrogen in the Solar 
System is the 14 N/ 15 N isotopic ratio. While good data exist for the Solar System, corresponding measurements in the 
interstellar medium are scarce and of low quality. 

Aims. Following the successful detection of the singly, doubly, and triply deuterated isotopologues of ammonia, we have 
searched for 15 NH_2D in dense cores, as a new tool for investigating the 14 N/ 15 N ratio in dense molecular gas. 
Methods. With the IRAM-30m telescope, we have obtained deep integrations of the ortho 15 NH2D (li,i — lo,i) li ne a t 
l/") ■ 86.4 GHz, simultaneously with the corresponding ortho NH2D line at 85.9 GHz. 

Results. The ortho 15 NH 2 D (li - 1 1) is detected in Barnard-lb, NGC1333-DCO+, and L1689N, while we obtained 
upper limits towards LDN1544 and NGC1333-IRAS4A, and a tentative detection towards L134N(S). The para line at 
£^ , 109 GHz remains undetected at the rms noise level achieved. The 14 N/ 15 N abundance ratio in 15 NH2D ranges between 

350 and 850, similar to the protosolar value of ~ 424, and likely higher than the terrestrial ratio of ~ 270. 

Key words. ISM clouds - molecules - individual object (Barnard-lb, L1689N, L134N(S), L1544, NGC1333-IRAS4A) - 
radio lines: ISM 
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RA 
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n(H 2 ) a 




(2000) 


(2000) 


( km s —1 ) 


(cm ~ 3 ) 


Barnard lb 


03:33:20.9 


31:07:34 


6.8 


3 x 10° 


NGC1333-IRAS4A 


03:29:10.5 


31:13:31 


7.2 


2 x 10 6 


NGC1333-DCO+ 


03:29:12.3 


31:13:25 


7.2 


1 x 10 6 


L1544 


05:04:16.6 


25:10:48 


7.4 


2 x 10 6 


L134N(S) 


15:54:08.6 


-02:52:10 


2.4 


2 x 10 6 


L1689N 


16:32:29.5 


-24:28:53 


3.8 


2 x 10 6 



1. Introduction Table 1. Source list 

h: 

5~ , Nitrogen chemistry is particularly interesting for un- 
derstanding the connection between the ISM and the 
formation of the solar nebula, because it is thought that 
the primitive atmospheres were nitrogen rich, as Titan 
remains today. Furthermore, the isotopic 15 N/ 14 N ratio 
has been measured in a variety of Solar System bodies, 
from the giant planets to the rocky planets, comets, and 
meteorites. The observed differences in nitrogen fraction- 
ation are used to understand how these bodies formed 

within the protosolar nebula. The combination of nitrogen a From Caselli et al. ( 2008) 
and hydrogen (D/H) isotopic ratios has been demonstrated 
to be a very effective way of understanding how the ice 

mantles wer e enr iched in deuterium and nitrogen. Aleon 01 the organic matter, enriched in lo N and deuterium, 
and Robert (|2004T ) have concluded that a fast condensation is needed in order to keep a significant fractionation in 

the solid material of the primitive Solar System. They 

Send offprint requests to: e-mail: maryvonne.gerin@ens.fr al SO have evaluated the exothermicity of the fractionation 

* Based on observations obtained with the IRAM 30 m tele- reactions for nitrogen to be 43 ± 10 K. The D fraction- 
scope. IRAM is supported by INSU/CNRS (France), MPG ation has not been inherited from the native prestellar 
(Germany), and IGN (Spain). core, but most likely occurred in the protosolar nebula 
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Table 2. Einstein coefficients, upper level energies and crit- 
ical densities for the range of temperatures considered in 
this work 




Fig. 1. Spectra of the — lo,i lines of o-NH 2 D (grey/red 
line) and o- 15 NH2D(black line) . The vertical scale is T m bin 
K, the horizontal scale is Vlsr m kms -1 . The o-NH 2 D 
spectra have been multiplied by 0.02 except for the L1689N 
spectra that has been scaled by 0.0125. 



(author?) (jRemusat et al. 20061 iGourier et al. 20081 ). yet 
the same physical and chemical processes are thought 
to operate in the prestellar cores and in the coldest 
regions of circumstellar disks. In the ISM, recent ob- 
servations show that, contrary to CO, nitrogen does 
not deplete from the gas phase in dense cores, except 
when the density rises significantly above 10 6 cm~ 3 . 
Nitrogen species can therefore be very significantly 
deuterated, with D/H fractio n ation of several tenths 
for N 9 p+ (iDaniel et al. 20071 ; iPagani et al. 20071 ) and 
NH 2 D (jCrapsi et al. 20071 ). Multiply de uterated ammonia 
in particular c an be very abund ant dGerin et al. 2 006; 
iLis et al. 2002a! iLis et al. 200l iRoueff et al. 2005t) . 
Nitrogen molecules will therefore be significant molecular 
reservoirs of deuterium. It is interesting to study whether 
they could also be enriched in 15 N, and whether signatures 
from an enrichment at an early evolutionary stage can be 
identified in primitive matter. 

High 15 N enhancements are measured both in HCN 
and CN cometa ry gases ()Bockelee-Morvan et al. 2008t 
2008 |). and in primitive carbonaceous me- 
15 N enhancements may have been present 



ISchulz et al 
teorites. High 



in the ammonia ices of the natal presolar cloud accord- 
ing to the fr actionation mechan ism proposed by Rodgers 
fe Cha rnley (|2008al : l2008bt I2004D and Charnley & Rodgers 
(2002). Nitrogen fractionation is not expected to be as effi- 
cient as deuterium fractionation in dense cores, yet signifi- 
cant departures from the elemental 14 N/ 15 N ratio may oc- 
cur in some molecules. As first shown by Terzievi a fe H erbst 
(2000), and develope d by Ch arnley & Rodgers (|2002f ) and 
Rodgers & Charnley (|2008al ), nitrogen fractionation in the 
gas phase may operate through ion-molecule reactions in- 
volving atomic or ionized nitrogen . Rodgers and Charnley 
(2008b{) have subsequently studied the possible role of 
neutral-neutral reactions involving 15 N and CN. Little ob- 
servational interstellar data are available. We have therefore 
started a survey of the main nitrogen-bearing interstellar 
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o-NH 2 D 


— lo,i 


85926.2780 


7.82e-6 


20.68 


4.2 10 b 


o- 15 NH 2 D 


li.i — lo,i 


86420.1959 


7.96e-6 


20.63 


4.2 10 6 


p- 15 NH 2 D 


— lo,i 


109284.9021 


1.61e-5 


21.18 


8.8 10 6 



species in 5 dense cores and a class source (Table [IJ. This 
paper reports the detection of o- 15 NH 2 D as the first result 
of this survey. 



2. Observations 

The microwave and far infrared spectra of 15 NH 2 D and 
15 NHD 2 have been recently investigated by Elkeurti et al. 
(|2008f ) and used to produce the corresponding line lists as 
supplementary data_j, while accurate line lists and partition 
functions for the 14 N isotopol ogues of the NH 3 family can 
be found in Coudert & Roueff (|2006l ) . These species are also 
independently listed in the Cologne Database for Molecular 
Spectroscopy (CDMS, Miiller et al. 2001, 2005), with small 
differences in the line frequencies due to different handling 
of the hamiltonians. 

We have chosen to search for the li^ — lo,i line of ortho 
15 NH 2 D since the corresponding NH 2 D line is very strong 
and both the sky transmission and telescope performances 
are excellent at 86 GHz. The frequency shift introduced 
by the 15 N substitution is small enough that the two iso- 
topologues can be observed with the same receiver tuning. 
The line frequencies (Elkeurti et al. 2008), Einstein A coeffi- 
cients, upper energy levels and critical densities are listed in 
Table [21 We have used the theoretical estimates of the criti- 
cal densiti es fro m the reduced mass ratio scaling of Machin 
& Roueff (|2006l ) for the NH 2 D-He values at 10 K, the tem- 
perature appropriate for the cold cores we have observed. 
However these values are probably too large when molecu- 
lar hydrogen is involved, as found in recent calculations of 
the NH3-H2 system by Valiron et al. (private communica- 
tion) . 

The observations have been performed with the IRAM- 
30m telescope, during three observing sessions in December 
2007, March 2008 and September 2008. We used the A100 
and B100 receivers in parallel, tuned to 86.2 GHz in order to 
detect o-NH 2 D and o- 15 NH 2 D with the same detector set- 
ting. The weather conditions were average, with 5 — 10 mm 
of water vapor (PWV) . The NH 2 D and 15 NH 2 D 
lines were observed simultaneously, with the J—l— 
lines of H 15 NC and H 13 CN at 86.055 GHz and 
86.338 GHz. We used the VESPA correlator, tuned 
to a spectral resolution of 40 kHz, and spectral bandpass 
of 40 MHz for each line. The data were taken using the 
wobbling secondary reflector, with a beam separation of 
240" . Telescope pointing was checked on nearby planets and 
bright radio quasars and was found accurate to ~ 3". Due 
to rather poor weather conditions during the September run 



1 Available at http : //library, osu. edu/ sites/msa/ suppmat/ 
v251 . il-2 .pp90-101/mmcl . txt 
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(high PWV and cloudy sky), the pointing accuracy was de- 
graded to ~ 5". Additional observations of the p- 15 Nti2D 
line at 109.3 GHz were obtained in March 2008. We only 
searched for this line towards Barnard-lb and detected no 
signal down to a rms noise level of 18 mK with 0.2 kms -1 
velocity resolution. For L134N(S), we combined the data 
with observations performed i n April 2005, as part of the 
dark cloud line survey project (Marcclin o et al. 20091) . The 
weather conditions were excellent (1-2 mm PWV) and the 
observations performed in the frequency switching mode. 

The data processing was done with the GILDAE0 soft- 
ware (e.g. Pety et al. 2005). We used the dec08b version of 
this software, which allows to correct for a minor bug in the 
frequency calibration during the observations. The IRAM- 
30m data are presented in main beam temperatures T m t>, 
using the forward and main beam efficiencies F c g and B c ff 
appropriate for 86 GHz, F c g—0.95 and -B o ff=0.78. The un- 
certainty in flux calibration is ~ 10%, as checked by the 
variation of the intensity of the strong 0-NH2D and H 13 CN 
lines in the spectrum. Linear baselines were subtracted. 

Because the nuclear spin of 15 N is 1/2, the 15 NH 2 D 
lines are split into fewer hyperfine components than NH2D 
which makes their detection more favorable. The hyper- 
fine structure of 15 NH2D is driven by the quadrupole mo- 
ment of the deuterium nucleus, which is much smaller than 
the corresponding value of 14 N. We have checked, by us- 
ing the nuclear q uadru pole constants provided in Garvey 
et al. (author?) (|l976h . that the resulting hyperfine split- 
ting is less than 50 kHz. We can thus safely assume that 
the spectrum reduces to a single component . As shown 
in Figure [U the 15 NH2D line is clearly detected towards 
Barnard-lb, and L1689N, while we obtained upper limits 
towards LDN1544 and NGC1333-IRAS4A and tentative 
detections towards NGC1333-DCO+ and L134N(S). 
The ratio of peak antenna temperatures of the NH2D and 
15 NH2D lines is 50 - 100, and the velocity agreement 
is excellent. Using the JPL and CDMS spectroscopy data 
bases, we have checked that no line of known interstellar 
molecules are expected within ±300 kHz from the 15 NH2D 
line frequency. The identification of the detected feature is 
therefore secure. 

The line parameters were estimated by fitting Gaussian 
profiles to the detected o- 15 NH2D lines. For o-NH 2 D we 
used the HFS routine implemented in CLASS, which al- 
lows to take into account the hyperfine components self- 
consistently. The opacity of the ortho NH 2 D line is moder- 
ate in all sources, with a total opacity for all lines ranging 
from - 1 to - 5 (Table©. 

3. Results 

3.1. NH 2 D and 15 NH 2 D 

Results of the fits and derived molecular column densities 
are listed in Table[3] As we are mostly interested in the ratio 
of column densities, we computed them under the simple 
assumption of a single excitation temperature. We used the 
excitation temperature derived from the NH2D fit for both 
isotopic species. The o-NH 2 D column densities are in good 
agreeme nt with previously published results for the same 
sources (|Roueff et al. 20051) . The [NH 2 D]/[ 15 NH 2 D] abun- 
dance ratio range from 360 to 810, with the largest value for 

2 See http : //www . iram . f r/IRAMFR/GILDAS 



3 

L1689N. This last source is an interaction region between a 
molecular outflow a nd a dense core, and as such may have 
peculiar properties (|Lis et al. 2002b| ). Given the error bars, 
the measured [NH2D]/[ 15 NH2D] ratio is comparable to the 
14 N/ 15 N protosolar ratio, as measured in Jupiter 
(450; Fouchet et al. 2004) and in osbornite-bearing 
calcium-aluminium-rich inclusions from meteorites 
(424; Meibom et al. 2007), and likely larger than 
the terrestrial abundance ratio (270). Although the 
uncertainty on the [NH2D]/[ 15 NH2D] ratio remains large, 
the cold prestellar cores L1689N and LDN1544 seem to have 
a larger ratio than Barnard-lb and NGC1333-DCO + . 

3.2. 15 N fractionation 

Nitrogen fractionation involves two main mechanisms in the 
gas phase: isotopic dependent photodissociation of molec- 
ular N2, principally at work in the atmosphere of Titan 
(Liang et al. 2007) and possible ion-molecule fractionation 
reactions occurring at low temperatures in cold dense cores 
as first measured by Adams & Smith (1981) and calcu- 
lated by Terzieva and Herbst (2000). In this latter case, 
involved endothermicities values range between a few K up 
to 36 K for exchange reactions involving 15 N, 15 N + , and 
15 NN. Selective photodissociation of N 2 and 14 N 15 N 
takes place at wavelengths between 80 and 100 nm, 
a range where cold dense cores are completely 
opaque. Then, this mechanism does not work in the 
present con text. Charley & Rodgers (120021 ) and Rodgers 
& Charnley (j2008al ) have investigated the nitrogen fraction- 
ation in their time dependent, coupled gas/solid chemical 
models. They conclude that 15 N-rich ammonia and deuter- 
ated ammonia can be frozen onto the ice mantles, provided 
all nitrogen is not converted into N2. The gas phase be- 
comes enriched at early times, before the complete freezing 
of the gas phase molecules. 

Additional fractionation reactions may be introduced 
such as those involving 15 N + with CN and NH3 and some 
neutral- neutral reactions between 15 N and CN (Rodgers 
& Charnley 2008b). However, none of these reactions has 
been studied in the laboratory and these schemes remain 
highly hypothetical. We have developed a gas phase chem- 
ical code, including ion- molecule fractionation reactions 
for carbon and nitrogen (jLanger 19921 : lLanger et al. 198~4l 
iTerzieva fc Herbst 2000h as well as a complete deuterium 
chemistry (Roueff et al. 2005). We have explicitly intro- 
duced D and 13 C on the one hand and D and 15 N on the 
other hand for NH„, HCN and HNC molecules, in order to 
directly compare the model results with the observations. 
The chemical network involves 302 chemical species and 
5270 reactions. The maximum number of carbon atoms in 
a molecule has been limited to 3. We have introduced the 
additional reaction channels arising from the inclusion of 
isotopic species. We have also preserved functional groups 
in dissociative recombination reactions such as : 

HCND+ + e -> HCN + D (1) 

HCND+ + e DNC + H (2) 

Note that the branching ratios of the dissociative re- 
combination of N 2 H + have been measured again by Molek 
et al. (2007) with the result that the channel towards N 2 
occurs with a probability of at least 90% 



4 Gerin et al.: Detection of interstellar NH2D 

Table 3. Line intensities and molecular column densities 
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T 
i 

mKkms" 
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iV 

10 11 cm" 2 


Barnardlb 


2.5 ±0.047 


0.79 


5.24 ±0.14 


6.0 ±0.5 


4.7 ±0.5 


42 ±9 


30 ±4 


0.67 


10 ±2.7 


47n +ivu 

4 ' U -100 


N1333-IRAS4A 


1.0 ±0.018 


1.38 


1.39 ±0.10 


5.0 ±0.5 


2.7 ±0.6 


±10 


< 30 




< 10 


> 270 


N1333-DCO+ 


1.3 ±0.015 


1.15 


1.71 ±0.05 


5.3 ±0.5 


2.4 ±0.4 


26 ±8 


14 ±3 


0.52 


6.7 ±2.5 


3601^ 


LDN1544 


2.3 ±0.016 


0.47 


7.05 ±0.05 


5.3 ±0.5 


4.1 ±0.5 


±7 


< 10 




< 5.2 


> 700 


L134N(S) 


2.2 ±0.033 


0.42 


4.75 ±0.10 


5.5 ±0.5 


2.4 ±0.4 


24 ±7 


10 ±2 


0.40 


4.5 ± 2 


530± 5 ^ 


L1689N 


5.3 ±0.030 


0.53 


6.98 ±0.02 


8.5 ±0.5 


3.4 ±0.5 


65 ± 17 


26 ±6 


0.37 


4.2 ± 1.5 


810l 6 J 



a a is the rms computed for the original spectral resolution of 40 kHz = 0.136 kms . 

b computed at LTE with the same T ex for 0-NH2D and o- 15 NH 2 D. T ex is derived from the HFS fit of the o-NH 2 D profile. 




Fig. 2. Prediction of the gas phase abundances relative to 
H2 (left) and 14 N/ 15 N abundance ratio (right) for the main 
nitrogen species. The models assumes a constant temper- 
ature of 10K, and increasing depletions with the gas den- 
sity, to mimic freezing out. The elemental abundance ratio 
14 N/ 15 N is set to 400. 

A calculation is shown in Figure [2] for typical dense core 
parameters, and assuming a 14 N/ 15 N abundance ratio of 
400, and an ionization rate of ( = 2 x lO -17 s _1 . The model 
predicts that the 15 N enrichment of ammonia is moderate 
in the gas phase, while a stronger enrichment is predicted 
for N2H" 1 ", and depletion for HCN and CN. Recent models 
by Rodgers and Charnley (2008a) obtain similar results for 
the gas phase abundances, the 15 N enrichment of ammonia 
being more efficient in the solid phase. 

4. Conclusions 

We report the detection of heavy deuterated ammo- 
nia, 15 NH2D, in three cold dense cores. The abun- 
dance ratio [NH 2 D]/[ 15 NH 2 D] is compatible with 
the 14 N/ 15 N protosolar value, and seems larger than 
the terrestrial value despite the remaining mea- 
surement uncertainties. While further observations are 
needed to improve the accuracy and test our chemical mod- 
els, ammonia and deuterated ammonia seem to be good 
probes of the 14 N/ 15 N ratio. Deuterated ammonia is par- 
ticularly interesting as it probes the coldest and densest 
regions of prestellar cores which are the reservoirs for the 
future formation of young stars and their associated proto- 
planetary disks. 
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